Search for variation of fundamental constants and violations of fundamental 

symmetries using isotope comparisons 
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Atomic microwave clocks based on hyperfine transitions, such as the caesium standard, tick with 
a frequency that is proportional to the magnetic moment of the nucleus. This magnetic moment 
varies strongly between isotopes of the same atom, while all atomic electron parameters remain 
the same. Therefore the comparison of two microwave clocks based on different isotopes of the 
same atom can be used to constrain variation of fundamental constants. In this paper we calculate 
the neutron and proton contributions to the nuclear magnetic moments, as well as their sensitivity 
to any potential quark mass variation, in a number of isotopes of experimental interest including 

001 1 OQ ft*7 ft ^ 

Hg and ' Rb, where experiments are underway. We also include a brief treatment of the 
dependence of the hyperfine transitions to variation in nuclear radius, which in turn is proportional 
to any change in quark mass. Our calculations of expectation-values of proton and neutron spin in 
nuclei are also needed to interpret measurements of violations of fundamental symmetries. 



I. INTRODUCTION 

Studies investigating possible space-time variation of 
the fundamental constants of particle physics have been 
given additional motivation in recent years by indications 
of variation of the fine-structure constant, a = e 2 /he, in 
quasar absorption systems 0-01 ■ Such a variation could 
be confirmed using complementary astrophysical [4] and 
laboratory |5j studies, but so far all laboratory studies 
have shown results consistent with zero variation (see, 
e.g. the reviews d, 0])- 

Theoretical models that attempt to unify the forces 
of the Standard Model (electroweak and strong nuclear) 
with gravity consistently predict that a variation in a 
would be related to a variation in the quantum chromo- 
dynamics (QCD) scale and the Higgs field. Such varia- 
tions could be measured as a change in the dimension- 
less ratio m q / Aqcd, where m q is the light quark mass 
(formally the average of the up and down quark masses, 
uiq = (m u + m,d)/2) and Aqcd is the quantum chromo- 
dynamics constant (which shows the position of the pole 
in the running strong coupling constant). The variation 
of m q / Aqcd has been estimated to be as great as 35 
times that of a, and consequently may be more suitable 
for investigation (see, e.g. |8l-[l2|). 

Measuring a change in m q j 'Aqcd directly is not pos- 
sible, so some effort must be made to calculate the re- 
lationship between variation in quark mass and that of 
an observable quantity. The dependence of proton and 
neutron magnetic moments on m q j Aqcd has previously 
been calculated [l3|, [l4[ . These calculations in turn can 
be used to obtain the sensitivity of nuclear magnetic mo- 
ments, fi, to variations in m q j Aqcd .15]. Nuclear mag- 
netic moments can be probed by the hyperfine structure 
of atomic clocks, which form the reference transition in 
clocks such as the caesium frequency standard. 

By comparing the hyperfine frequencies of isotopes of 
the same atom, a stringent limit can be placed on vari- 
ations in m q j Aqcd- The comparison is insensitive to 
a-variation since this dependence comes only from the 



electronic structure, which is common to both. For the 
same reason, difference measurements between isotopes 
of the same ion have strongly reduced sensitivity to pres- 
sure shifts, second-order Zeeman shifts, and second-order 
Doppler shifts fl6i | . Experiments currently underway in- 
clude comparison of 201 H g an d 199 Hg clocks 0, E3' an d 
87 Rb versus 85 Rb clocks [18| ; for limits from other mi- 
crowave clock experiments see the review [l9j . Other 
candidate isotopes are also considered in this paper. In 
this work we assume that Aqcd is constant, calculating 
the dependence on the small parameter m q . This is not 
an approximation: rather it means that we measure m q 
in units of Aqcd- Therefore Sm q /m q should be under- 
stood as a relative change in m q j Aqcd- 

In this work we have two aims. First, we find expec- 
tation values of proton and neutron spin for isotopes of 
nuclei of experimental interest. These are needed in order 
to calculate nuclear electric dipole moments (EDMs) 



(1) 



due to proton and neutron EDMs {d p and d n , respec- 
tively). Here / is the nuclear spin and s p = s n = 1/2. 
Similarly we may calculate constants of the parity non- 
conserving (PNC) nuclear-spin-dependent interaction be- 
tween electrons and the nucleus (see, e.g. [13, Hll) using 



(2) 



where (s) = (s) z for the maximal projection I z = I; we 
omit the subscript z for brevity. Values of (s p ) and (s n ) 
are presented for many nuclei of experimental interest in 
Tables U and [TTJ In these tables the "BP column rep- 
resents our preferred value, while the difference between 
BI and BII gives an indication of the uncertainty in the 
calculation. 

Our second aim is to use these expectation values to 
calculate the sensitivity of the nuclear magnetic moment 
to variation of quark mass, denoted k. We then improve 
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our estimates of k by including the variation in proton 
and neutron-spin expectation values induced by quark 
mass variation, as well as nuclear radius effects. Results 
for nuclei of experimental interest are presented in Ta- 
ble IIIII We identify potential candidate nuclei for exper- 
imental study where there is high sensitivity to quark- 
mass variation and good stability of calculations across 
the various theoretical models. 



II. NUCLEAR MAGNETIC MOMENT 

The first step in our calculations is to disentangle the 
proton and neutron contributions to the nuclear mag- 
netic moment (in units of the nuclear magneton (J-n = 



/i = g p (s p ) +g n (s n ) + (l p ) , 



(3) 



where {s p } and (s n ) are the expectation values of the total 
proton and neutron spins, while (l p ) is the orbital angu- 
lar momentum of an unpaired valence proton (if there is 
one). In this paper we treat the case where the nucleus 
has either one valence neutron or one valence proton: 
because of nuclear pairing even-even nuclei have no spin 
(and hence no magnetic moment), while odd-odd nuclei 
are often unstable (of course with many exceptions) and 
tend not to be used in microwave clocks. We begin by 
considering the contribution of the valence nucleon alone. 



Model A 

Assuming all other nucleons in the nucleus are paired 
and ignoring polarization of the core (which will be con- 
sidered in model B), the spin /, and hence magnetic mo- 
ment of the nucleus, is entirely due to the total angular 
momentum of that external nucleon: I = j = I + s. One 
obtains the standard Schmidt formula for the magnetic 
moment: 



fi° = g s (s z ) +gi (l z )° 

if j=l + \ 

2(i+l) ' J 2 



(0° = { 1 



(4) 
(5) 
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The gyromagnetic factors are gi = 1, g s = g p = 5.586 
for a valence proton and gi = 0, g s — g n — —3.826 for a 
valence neutron. It is these Schmidt values for (s p ) and 
(s n ) that are presented as the model A results in Tables U 
andHH 

Because light-quark masses are small, the g-factors are 
sensitive to m q / Aqcd mainly via 7r-meson loop correc- 
tions to the nuclear magnetic moments, which are pro- 
portional to 7r- meson mass ~ \Jfn q ^-QCD [HI- Be- 
cause m n = 140 MeV, the quark-mass contribution can 



be significant even in the chiral limit m q — 0. Full details 
are given in [l3j (see also [3); the final results are 



8g P 

9n 



-0.087^ 



0.013- 



8m , 



-0.118^+0.0013 



m s 
5m s 
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where variations in m q and m s (strange quark mass) are 
considered separately. Note that the g-factors are far 
more sensitive to m q than m s . If one assumes the quark- 
masses as proportional to the vacuum expectation value 
of the Higgs field, then Sm q /m q should be approximately 
equal to Sm s /m s . The total variation Sg/g would then 
be the sum of the variations in m q and m s . 

In is convenient at this stage to define parameters K p 
and K n by 



5n 
I* 



Sgp 
9 P 



Kn 



9n 



where, according to ©, 



K p = 



9p {s P } 



and Kr, 



9n (s n ) 



(7) 



(8) 



For example, for a valence proton, model A gives K p — 

9p ( s p)° /m° anc l Kn — 0. Eq. (J?]) shows the connec- 
tion between the spin expectation value and sensitivity 
to variation of fundamental constants. 

Combining Eqs. ([5]) and (J7J one may obtain the sen- 
sitivity of the magnetic moment /i to variation in light 
and strange quark masses for any nucleus with a single 
valence nucleon (more details may be found in [15[). The 
result is characterized by the sensitivity coefficients K q 
and k s , defined by 



5[i 



Sm a 



+ Ks 



8m« 



(9) 



where 



K q = -0.118 K n - 0.087 K p 
k s = 0.0013 K n - 0.013 K p . 

In this paper we are interested in finding large differences 
in k = K q + k s between different isotopes of the same 
nucleus. With this in mind, we present in Table HTT1 the 
k obtained using model A, via Eq. ([5]), for all isotopes 
of atoms commonly used in microwave clocks. We defer 
discussion of the results to Section HVl 



Model B 

The shell model of the previous section is known to 
overestimate the magnetic moment determined from ex- 
periment. This can be understood as polarisation of the 
non-valence nucleons which acts to reduce the magnetic 
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TABLE I. Average values of (s p ) and (s„) for heavy nuclei 
with a valence proton. A: Schmidt values (no neutron contri- 
bution); BI: our preferred fit of experimental magnetic mo- 
ment assuming separate conservation of j p , j n , and total I and 
s; BII: fit \i assuming minimal transfer of spin from protons 
to neutrons. 







/a \ 
\ S P> 








Isotope 


A 


BI 


BII 


BI 


BII 


133 Cs 


—0.389 


—0.286 


—0.297 


—0.103 


—0.092 


135 Cs 


—0.389 


—0.268 


—0.281 


—0.121 


—0.108 


137 Cs 


—0.389 


—0.255 


—0.269 


—0.134 


—0.120 


131 Cs 


0.500 


0.351 


0.367 


0.149 


0.133 


129 Cs 


0.500 


0.345 


0.362 


0.155 


0.138 


127 Cs 


0.500 


0.341 


0.358 


0.159 


0.142 


125 Cs 


0.500 


0.335 


0.353 


0.165 


0.147 


139 Cs 


—0.389 


—0.272 


—0.285 


—0.116 


—0.104 


123 Cs 


0.500 


0.332 


0.35 


0.168 


0.150 


i39 La 


—0.389 


—0.262 


—0.276 


—0.127 


—0.113 


137 La 


—0.389 


—0.273 


—0.285 


—0.116 


—0.104 


'2'2-i Fj . 


0.500 


0.188 


0.221 


0.312 


0.279 


221p r 


-0.357 


-0.272 


-0.281 


-0.085 


-0.076 


225 Fr 


0.500 


0.176 


0.211 


0.324 


0.289 


227p r 


0.500 


0.346 


0.363 


0.154 


0.137 


203 rpj 


0.500 


0.361 


0.376 


0.139 


0.124 


205 rpj 


0.500 


0.363 


0.377 


0.137 


0.123 


201 rpj 


0.500 


0.359 


0.374 


0.141 


0.126 


199 rpj 


0.500 


0.358 


0.373 


0.142 


0.127 


197rpj 


0.500 


0.356 


0.371 


0.144 


0.129 


195 rpj 


0.500 


0.356 


0.371 


0.144 


0.129 


193 rpj 


0.500 


0.357 


0.372 


0.143 


0.128 


207rpj 


0.500 


0.391 


0.403 


0.109 


0.097 



moment. Using the experimental value of the magnetic 
moment we can estimate this core polarisation, and thus 
improve our estimates for (s p ) and (s n ). 

There are many ways to enact a reduction in magnetic 
moment from the Schmidt value /i° to the experimen- 
tal value /z. The most efficient means is to assume that 
the spin-spin interaction transfers spin from the valence 
proton (neutron) to core neutrons (protons): 

((Sp) ~ (S P )°) = -««„> - <*n>°) = (10) 

g P - g n 

where (s n )° and (s p )° are the Schmidt values from model 
A (one will necessarily be zero). The denominator 
9p~ 9n = 9.412 is a large number, so the required change 
in (s p ) and (s n ) to obtain the experimental /j, is mini- 
mal. The transferred spin goes from contributing with a 
gyromagnetic factor of g p = +5.586 to g n — —3.826. 

Of course nature may proceed in an entirely different 
fashion, for example by transferring spin of a valence pro- 
ton to its orbital angular momentum, or in a rather more 
unlikely fashion by transferring valence proton spin to 
neutron orbital momentum. Even under extreme scenar- 
ios, it was shown in [15| that the difference in the correc- 
tion to the m g -dependence varied by less than 10%. This 
stability of the final result to the assumptions enables us 



TABLE II. Average values of (s p ) and (s n ) for nuclei with 
a valence neutron. A: Schmidt values (no proton contribu- 
tion); BI: our preferred fit of experimental magnetic moment 
assuming separate conservation of j p , j n , and total I and s; 
BII: fit /i assuming minimal transfer of spin from protons to 
neutrons. 
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A 


B I 


B II 


B I 


B II 


129 Xe 

Ac 


u.ouu 




n ^7Q 


U. LOO 


n 1 91 


131 Xe 


-0.300 


-0.246 


-0.252 


—0.054 


-0.048 


127 Xe 


0.500 


0.332 


0.35 


0.168 


0.150 


133 Xe 


—0.300 


—0.26 


—0.264 


—0.04 


—0.036 


135 Xe 


—0.300 


—0.271 


—0.274 


—0.029 


—0.026 


iYi Yh 


n i R7 


n im 

— U. lOU 


— U.1J1 


n ni 7 


n ni ^ 


173 Yb 




— n 1 1 r 


— n i a^k 

o. l^O 




—0.214 


169 Yb 




— n i X7 


0. 1UO 


—0.252 


—0.226 


175 Yb 


n ^nn 


n i si 


n 91 ^ 


n ^i q 


n 98^ 


167 Yb 


U .OO 1 


— n 9RQ 


— n 97R 


— n ORR 
u.uoo 


— n n7Q 


163 Yb 


0.500 


0.317 


0.336 


0.183 


0.164 


165 Yb 


—0.357 


—0.252 


—0.263 


—0.106 


—0.094 


161 Yb 


0.500 


0.311 


0.331 


0.189 


0.169 


199 H> 


— n 1 fi7 


— n 1 'ii 


— n 1 'i^ 

O. ItJO 


— n ni fi 


— n ni a 


201 Hg 


0.500 


0.339 


0.356 


0.161 


0.144 


203 H 


—0.357 


—0.296 


—0.302 


—0.062 


—0.055 


197 He 


—0.167 


—0.154 


—0.155 


—0.013 


—0.012 


195 He 


—0.167 


—0.155 


—0.156 


—0.011 


—0.010 


193 Hg 


0.500 


0.347 


0.363 


0.153 


0.137 


191 Hg 


0.500 


0.346 


0.362 


0.154 


0.138 


189 Hg 


0.500 


0.345 


0.361 


0.155 


0.139 


205 Hg 


-0.167 


-0.162 


-0.163 


-0.004 


-0.004 


225 Ra 


0.500 


0.360 


0.375 


0.140 


0.125 


223 Ra 


-0.300 


-0.196 


-0.207 


-0.104 


-0.093 


227 Ra 


-0.300 


-0.116 


-0.135 


-0.184 


-0.165 


229 Ra 


0.500 


0.213 


0.243 


0.287 


0.257 


211 Rn 


-0.167 


-0.162 


-0.163 


-0.004 


-0.004 


209 Rn 


-0.357 


-0.294 


-0.301 


-0.063 


-0.056 


207 Rn 


-0.357 


-0.292 


-0.299 


-0.065 


-0.058 


205 Rn 


-0.357 


-0.290 


-0.297 


-0.067 


-0.060 



to ignore the detailed nuclear forces and use "heuristic" 
fitting of \i to obtain values of (s p ) and (s n ). 

The preferred method of fitting /u, in [ll| is to assume 
that the total z-projection of proton and neutron angular 
momenta, j pz and j nz , are separately conserved, and that 
total spin and orbital angular momenta z-projections, 
(s p ) + (s n ) and (l p ) + (l n ), are also separately conserved 
(corresponding to neglect of the spin-orbit interaction). 
The final result gives values somewhere in between those 
of the "extreme" assumptions. Then 

(s z }° = (s P ) + (s„) 
(jpz) = (l P ) + (s P ) 

where (jpz) — I for valence proton and is zero for a va- 
lence neutron and (s z )° is the Schmidt value for the spin 
of the valence nucleon ([5]) . Manipulation of these expres- 
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sions with ([3]) gives 



For a valence proton: 



M~ Op*) ~ (flp ~ !) ( s z) 



(S n ) = 

g n - g P 

(Sp) = (s z )° - (s n ) . 



(11) 



These neutron and proton spins are presented for isotopes 
of elements used (or proposed) in microwave clocks as 
the BI column of Tables HI and ITT1 The "minimal change" 
model (|10p is presented as the BII column, and serves 
to indicate the expected uncertainty in our calculation. 
Sensitivity coefficients k calculated using values of (s p ) 
and (s n ) given by (|TT|) (corresponding to model BI) are 
given in the model B column of Table Mil 



K K = 



(fln ~9p + 1) ( s n) 



I* 



while for a valence neutron: 

(fJ P ~9n ~ 1) (s P ) 



K bn = 



(15) 



(16) 



The hnal sensitivity of our magnetic moment on m q is 
given in model C by 

k= -0.117 K n - 0.100 K p - 0.11 K b . (17) 

Calculations of k according to model C are presented in 
Table ITTT1 along with results of the other two methods. 



Model C 



III. HYPERFINE TRANSITIONS 



To calculate the expectation values (s p ) and (s n ), 
model B is the best we can do without venturing into 
highly model-dependent and complex nuclear calcula- 
tions. However, we should also estimate how these spin 
expectation values will change under a variation of quark 
mass, and include this in our calculation of K q . Under the 
assumption of conservation of nuclear spin (model B), we 
may define a parameter b in the case of a valence proton 
by 



(s„) = b (s z )° 

(s p ) = (l-b) (a,) 



(12) 



The coefficient b is determined by the spin-spin interac- 
tion and is estimated in model B by (fTTj) . In the case of 
a valence neutron, we should define the small parameter 
6 by (s„) = (1 - b) (s z ) r 

Under a variation of m 
according to 



instead of (fT2|) . 

q , the parameter b will change 



Sb = _ Qn Sm 1 
b m a 



(13) 



Full details are presented in [15!]. Briefly, the spin-orbit 
splitting remains finite in the chiral limit, and so it only 
has a weak dependence on m q / Aqcd- The spin-spin in- 
teraction is calculated in perturbation theory (to leading 
order) under the assumption that the major dependence 
comes from %- meson exchange (1/3 of the spin-spin in- 
teraction) and p- meson exchange (the remaining 2/3). 
If one instead assumes that the spin-spin interaction is 
dominated entirely by 7r-meson exchange, the coefficient 
in (Tr3"|) . —0.17, is not too different. 

We include the effect of variation of spin expectation 
values on magnetic moment by generalising our previous 
formula (0 to obtain 



To date, the most accurate probe of variation of fun- 
damental constants in the laboratory comes from atomic 
clocks. Microwave clocks based on the hyperfine interac- 
tion, such as the caesium fountain from which the SI time 
unit is defined, have achieved precision at parts in 10 15 
at several laboratories worldwide. To find the sensitivity 
of such clocks to potential variation of fundamental con- 
stants, we start with the Hamiltonian of the hyperfine 
transition which can be expressed as 



H, 



Ms 



AI J 



(18) 



where J and J are the nuclear and electron angular mo- 
menta, respectively. The A coefficient can be expressed 
as 



A = const 



h 2 



[a 2 F rel (Za) (1 - S h )] ^— (19) 



where Shifn, Za) is the small effect of finite nuclear ra- 
dius r„ [23| and F re \(Za) is the effect of relativistic 
corrections to the electron wavefunction at the nucleus. 
When comparing hyperfine transitions of two different 
nuclei, sensitivity to variation of fundamental constants 
will arise due to differences in fi, F IC \, and 5^. 

F Te i and Sh are sensitive to details of the electron wave- 
function, but are almost identical for different isotopes of 
the same element. Nevertheless, when comparing hyper- 
fine transitions across different elements F rc i and Sh may 
be included along with p,. Sh and p are mostly sensitive 
to quark mass, while -F rc i is sensitive to variation in the 
value of a. Therefore the effect of variation of these fun- 
damental constants on a hyperfine transition frequency 
u> (when measured against another hyperfine frequency) 
can be expressed 



Slu Sol . Sm„ 

— = A ro i h [K q + k hq ) 

u> a m„ 



(20) 



Here K [0 \ and kh q represent sensitivities of w to F Te \ and 
Sh, respectively. 
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TABLE III. Summary of magnetic moment sensitivity to quark- mass variation, k = K q + k s . Columns A, B and C refer to 
the different methods as described in Section II: model C is the best value, while the difference between B and C gives an 
indication of the theoretical uncertainty. The difference in k values for different isotopes of a particular element (relative to 
the first isotope) are given for models B and C. The isotope entries are given in order of half-life length, beginning with stable 
isotopes. 
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n nnft 

—u.uyo 


OS 


7 /9+ 


+Z.00Z 


otaDie 


n 1 97 

U. 1Z / 


n ri/i a 

U.U44 


n nn7 






135 n a 

OS 


7 /o+ 
7/2 T 


1 O 7QO 

+ Z. (OZ 


3 x 1U y 


U.1ZY 


n no r. 
U.Uoo 


n nnc 
— U.UOO 


n nnn 
— U.UU9 


n m o 
— U.Ulo 


Os 


7/2 T 


+ Z.001 


on i 7 , T 
3U.1 < y 


U.lz/ 


n non 

u.uzy 


n ni ^ 
— U.U14 


n ni c 
— U.Ulo 


n noo 
— U.UZZ 


131 r*„ 

OS 


o/z 


+o.oo(zj 


9.09 a 


— U.UOO 


n no7 
— U.U3 I 


n nno 
U.UUZ 


n no i 
— U.U81 


n nnr. 
—U.UUo 


129 rio 

OS 


1/2 T 


+i.4yi 


oo o u 
oz.o n 


— U.1UU 


n noo 
—U.Uoo 


n ni o 
U.Ulo 


n i 07 
— U.lz/ 


n nnft 
U.UUo 


Os 




i i a cn 
+ 1.459 


col, 

o.z n 


— U.1UU 


n noo 
— U.Uoz 


n m o 
U.Ulo 


n i oc 
— U.lzo 


n m i 
U.U11 


125 rio 
Os 


I/ 2 


i i /) no 
+ 1.4U9 


45 m 


n i nn 
— U.1UU 


n no i 
— U.Uol 


n noT 
U.Uz 


n i oc 
—U.lzo 


n non 
U.UzU 


139 r 
v^S 


7 /9+ 


-L9 fiQfi 
-rZ.oyo 


Q A m 


1 97 
U. 1Z / 


0. 037 


— U.UUO 


— U.UU / 


— U.U1U 


123 p„ 
OS 


1 /9+ 
1/2 


+ l.o / / 


O.o m 


n i nn 
— U.1UU 


n n©n 
— U.UoU 


n HQQ 

U.Uoo 


n 1 9/i 

— U. 1Z4 


n noft 
U.UZO 


La 


7 /0 + 

7/2 T 


1 O 700 

+z. /oo 


stable 


n i 07 
0.12/ 


n noo 
U.Uoz 


n m n 
— U.U1U 






La 


7/2 


i o cntr 
+ Z.090 


X 1U y 


n i 07 
0.1Z( 


n no7 
U.Uo ( 


n nno 
— U.UUo 


n nnc 
U.UUo 


n nri7 
U.UUf 


Y b 


1/2 


+0.494 


stable 


—0.117 


—0.116 


—0.084 






173 vu 

Y b 


c /o — 

0/2 


—0.648 


stable 


—0.117 


—0.125 


—0.467 


—0.009 


—0.383 


Y b 


7/2^ 


—0.635 


32.0 d 


—0.117 


—0.126 


—0.494 


—0.010 


—0.410 


1 75 \n a 

Y b 


7/2" 


+0.768 


4.18 d 


—0.117 


—0.126 


—0.510 


—0.010 


—0.426 


167 Y - 


5/2~ 


+0.623 


17.5 m 


-0.117 


-0.113 


0.018 


0.003 


0.102 


165 Yb 


5/2 _ 


+0.478 


9.9 m 


-0.117 


-0.112 


0.093 


0.004 


0.177 


161 Yb 


3/2 _ 


-0.327 


4.2 m 


-0.117 


-0.103 


0.430 


0.013 


0.514 


iyy Hg 


1/2" 


+0.506 


Stable 


-0.117 


-0.116 


-0.087 






201 Hg 


3/2- 


-0.560 


Stable 


-0.117 


-0.110 


0.156 


0.006 


0.243 


203 Hg 


5/2- 


+0.849 


46.8 d 


-0.117 


-0.115 


-0.048 


0.001 


0.039 


197 Hg 


1/2- 


+0.527 


64.1 h 


-0.117 


-0.116 


-0.093 





-0.006 


195 Hg 


1/2" 


+0.541 


9.9 h 


-0.117 


-0.116 


-0.097 





-0.009 


193 Hg 


3/2- 


-0.628 


3.80 h 


-0.117 


-0.111 


0.114 


0.005 


0.202 


191 Hg 


3/2- 


-0.618 


49 m 


-0.117 


-0.111 


0.120 


0.005 


0.207 


iB9 Hg 


3/2- 


-0.609 


7.6 m 


-0.117 


-0.111 


0.125 


0.005 


0.212 



The sign of the magnetic moment is unknown: we have assumed positive sign. 



For hyperfine transitions of s-wave electrons the fol- 
lowing formulae may be used: 



^hq 



-0.3 x 



(2 7 - l)S h 
S h 



K Ie 



(Za) 2 (12 7 2 - 1) 
7 2(4 7 2 _ !) 



where 7 = yfl — (Za) 2 and 
(21) 5 h « 1.995 (Z^/dBf- 1 



(22) 



(23) 
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TABLE IV. Sensitivity to quark mass variation due to nuclear 
size contribution, kh q , and sensitivity to a variation, K le \, 
for atoms and ions with s-wave valence electrons (these are 
defined by Eq. (Sty). 



Atom 




-Krcl 


Rb 


-0.003 


0.34 


Cd+ 


-0.005 


0.6 


Cs 


-0.007 


0.83 


Yb+ 


-0.014 


1.5 


Hg+ 


-0.023 


2.28 



Comparison with numerical calculations shows that 
Eq. ([2"Tj) tends to underestimate the a-dependence by a 
small amount (less than ~ 15%), while (|23l) includes a fit- 
ting parameter to the calculations that makes it accurate 
to within a few percent for all tested ions. In Table ITVl we 
present the results of previous numerical calculations of 
K re \ [l5[ and kh q [23j]. In the case where the microwave 
clock is based on hyperfine splitting of a p-wavc electron, 
such as proposed Al and Ga clocks, both corrections may 
be neglected. 



The same argument suggests caution in interpretation 
of the 201 Hg - 199 Hg comparison [l6j]. The sensitivity 
factor of K q = 0.24 predicted by model C represents a 
large departure from the model B value of 0.006. A more 
satisfactory result from this perspective is obtained from 
the 87 Rb - 85 Rb clock comparison, where models B and 
C give k = —0.16 and —0.08, respectively. The varia- 
tion between the models is not so strong, yet the final 
sensitivity factor, k — —0.08, is still reasonably large. 

We have also included a number of comparisons be- 
tween isotopes of different atoms in Table [Vj 5 k defines 
the sensitivity of a ratio, lo\/u2, of hyperfine transition 
frequencies of two isotopes to variation of quark mass: 

, . 5vn q 5m q 

- — - — r = (Ki-K 2 ) — -=5k — -. (24) 
(wi/w 2 ) m q m q 

These comparisons exhibit strong sensitivity as well as 
stability across the models. Once again, model C repre- 
sents our most complete calculation, while the difference 
between B and C gives an indication of the theoretical 
uncertainty. The final result includes the contribution of 
nuclear radius, 5k = SC + 5kh q . 



IV. DISCUSSION 

Using Table IIIII we may identify good isotopes for mi- 
crowave clock studies of variation of fundamental con- 
stants. The largest difference within a single element, 
SC in the table, is seen in ytterbium. Comparison of 
clocks based on hyperfine transitions in 161 Yb and 169 Yb 
should yield a relative sensitivity of 5k = 0.924. How- 
ever one should be cautious here. Because the valence 
nucleon is a neutron, the lowest order result (model A) 
is zero. The entire effect comes from the contribution of 
polarisation of the core protons (model B) and changes of 
this polarisation due to quark mass variation (model C). 
Thus, despite the apparent stability of (| 13[) as described 
earlier, the strong change in sensitivity k between models 
B and C should be treated with caution. 



TABLE V. Comparisons of sensitivity to quark-mass varia- 
tion, k, for isotopes of different atoms (see Eq. ([24])). Our 
final result is 5k = SC + 5kh q - 



Comparison: 


SB 


SC 


Skhq 


5 k 


133 Cs _ 199 Hg 


0.160 


0.095 


0.016 


0.111 


85 Rb - 199 Hg 


0.220 


0.151 


0.020 


0.171 


133 Cg _ 173 Yb 


0.169 


0.474 


0.007 


0.481 


133 Cs _ 169 Yb 


0.170 


0.501 


0.007 


0.508 


85 Rb - 173 Yb 


0.229 


0.531 


0.011 


0.542 
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